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CLAIMS 

1 . Method for estimating the movement between two numerical images h and b 
with luminance Y1 and Y 2 for generating for each point of coordinates x,y of the 

5 image I2 a movement vector d(x,y)=(d Xl d y ) so as to form an image I2 from the 
image h with luminance Y 2 ( x >y) = Yi(x-d x ,y-dy) which is an approximation of the 
image l 2l characterised in that it comprises the following steps : 

(a) defining an initial model of finished elements comprising a meshing 
whose nodes are points of the image I2, a movement vector at each node of said 

10 meshing, and an interpolation formula for calculating the value of the movement 
vector of each point of the image l 2 from the values of the movement vectors of 
the nodes of the mesh to which it belongs, 

(b) globally optimising the values of all the movement vectors of the model 
according to a differential method, 

15 (c) calculating a variation E between the image I2 and the image l 2 for each 

finished element or mesh, 

(d) carrying out a finer meshing on a discrete fraction of all the finished 

elements determined according to a criterion relating to the variations E and 

allocating a movement vector to each new meshing node, 
20 (e) repeating the steps (b), ( c ) and (d) on the model of finished elements 

obtained at the end of the preceding step (d) until a stoppage criterion is satisfied. 

2. Method according to claim 1, characterised in that, so as to carry out a finer 
meshing on a discrete fraction of all the finished elements in step (d), said set of 

25 finished elements is classified in the decreasing order of their variations E and the 
X first finished elements of this classification are subdivided into smaller finished 
elements, X representing a predetermined fraction of the number of finished 
elements of the set. 

30 3. Method according to claim 1 , characterised in that, so as to carry out a finer 
meshing on a discrete fraction of the set of finished elements in step (d), the set of 
variations E calculated in step (c) is compared with a threshold variation which 
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depends on the size of the finished element in question, and the finished 
elements whose variations E are greater than the threshold variation are 
subdivided into smaller finished elements. 

4. Method according to «=*e*e$ claim^ 1 characterised in that said stoppage 
criterion is a predetermined number of finished elements constituting the model of 
finished elements to be reached at the end of step (d). 

5. Method according to ©oe-ef-claim£ 1 te-^ characterised in that said stoppage 
criterion of step (e) is satisfied when the variations E of the set of finished 
elements of the model obtained at the end of the preceding step (d) are smaller 
than a functional threshold variation which depends on the size of the finished 
elements in question. 

6. Method according to ©Re-ef claim^ 1 te-§r characterised in that in addition for 
each numerical image U and fe, a set of R images l[ with a level of resolution r 
and luminance with r taking the values (0.....R-1) and i the values 1 and 2 is 
defined, the images I? and l° corresponding to the numerical images U and l 2 , 
and in that the steps (b) to (e) are carried out for each resolution level r from the 
level r=R-1 to the level r=0. 

7. Method according to claim 6, characterised in that the sets of R images with 
resolution level r are obtained by filtering the images h and l 2 along the two 
directions x and y using a low-pass filter with a pulse response h„ , each image l[ 
being defined by the following equation : 
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with M a natural integer. 
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8. Method according to claim 7, characterised in that the pulse response h r n is 
defined as follows : 
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s„ = 2B.sinc(27tB r n) = 2B 



sin27tB r n 
27rB r n 



B being a natural integer. 



9. Method according tc^efre~Of*the-pp 
initial movement vectors are nil vectors. 




ftng-etem^ characterised in that the 




characterised in that the 



variation E between the image l 2 and the image l 2 for each finished element e is 
defined by the following equation : 



1 1 . Method according to 0flo^Mh O" pr ooo d in g " Otarms > characterised in that the 
interpolation formula for calculating the value of the movement vector of a point P 
of coordinates (x,y) in the image l 2 belonging to the finished element e with 
vertices Pj, Pj and P k with respective coordinates (Xj.yO, (Xj.yj) et (x k ,yk) is the 
following : 



e= x DF ° 2 ( x .y) 



(x.y)ee 



where DFD(x.y) = Y 2 (x,y) - Yi (x-d x ,y-d y ) 




3(x,y)= x^ e ( x .y)*i.y.) 



I=i.j.k 

where is a function of the form 



(x,y) = a,+p,x + Yiy 

'Z^. ( x .y) = 1 



(x,y)ee 



%(x,y) = 0 



(x,y)ee 
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12. Method according to efte-ef- claimgf 1 to 14 ^ characterised in that the 
differential method for optimising the movement vectors is the Gauss-Newton 
method. 

5 

13. Method according to ©fte-ef claimp 1 « t o 4 4 » characterised in that the 
differential method for optimising the movement vectors is the Marquardt 
extension of the Gauss-Newton method. 

"io 14. Method according t(^ofte-e f""thoT3TO ^ characterised in that a 

compactness constraint is imposed on each finished element at the time of 
optimising the movement vectors of the model of finished elements, said 
constraint consisting of preventing the compactness of each finished element 
from tending to zero. 
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15. Method according to claim 14, characterised in that the compactness 
constraint on a finished element e with vertices Pj, Pj, P k and compactness 
C (PhPj.Pk) is defined by the following equation : 



20 where c^, c^, et represent the movement vectors of the vertices Pj, 

Pj, Pk during the optimisation step, and K is a compactness parameter. 

16. Method according to claim 14 or 1- 6? characterised in that the optimisation of 
the movement vectors under constraints on the finished elements is resolved by 

25 the increased Lagrangian technique. 

17. Method according to claim 16, characterised in that the constraints are used in 
a linearised form in the increased Lagrangian technique. 



30 18. Method according to efte-cff-clainrtf'^ 13 oi T7, characterised in that the 
methods for optimising the movement vectors use an LDL profile technique. 
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19. Application of the method for estimating movement between two numerical 
images according toeae^frclaimgf 1 4e-48-for coding images, characterised in that 
the fractional subdivision of the meshing carried out in step d) of the movement 
estimation method is associated with a partially quaternary tree in which each 
level represents a meshing level and each node represents a triangle of the given 
level, and in that what is generated is a binary train describing said tree. 

20. Application of the method for estimating movement between two numerical 
images according to claim 19, characterised in that the movement vectors 
associated with each node of said tree are encoded differentially with respect to 
the movement vectors of their father node when the latter exists and are ordered 
in said binary train along a width passage of said tree. 

21. Application of the method for estimating movement between two numerical 
images according to ©fte-ef claimjf l-te-^for decoding images, characterised in 
that the fractional subdivision of the meshing carried out in step d) of the 
movement estimation method is associated with a partially quaternary tree in 
which each level represents a meshing level and each node represents a triangle 
of the given level, and in that said tree is generated from a binary train of encoded 
data describing said tree. 

22. Application of the method for estimating movement between two numerical 
images according to claim 21 , characterised in that the encoded data relating to a 
given level of the tree are collectively regrouped in the binary train so as to 
generate the tree level by level as the train is read. 

23. Application of the method according to ®fte-e^foeT5recediT^^ to at least 
one of the ranges belonging to the group of the following ranges : 

- compression of sequences of images, and 

- compression of data in spaces larger than 2. 



